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Abstract. In this article we study a system of equations that is known to ex- 
tend Navier-Stokes dynamics in a weU-posed manner to velocity fields that are 
not necessarily divergence-free. Our aim is to contribute to an understanding 
of the role of divergence and pressure in developing energy estimates capable of 
controlling the nonlinear terms. We address questions of global existence and 
stability in bounded domains with no-slip boundary conditions. Even in two 
space dimensions, global existence is open in general, and remains so, primarily 
due to the lack of a self-contained energy estimate. However, through use 
of new coercivity estimates for the linear equations, we establish a number 
of global existence and stability results, including results for small divergence 
and a time-discrete scheme. We also prove global existence in 2D for any initial 
data, provided sufficient divergence damping is included. 



1. Introduction 

The zero-divergence constraint and the associated pressure field are the source 
of both difficulties and benefits in the study of the Navier-Stokes equations for 
the flow of viscous incompressible fluids. On one hand, the divergence constraint 
complicates analysis and approximation in a number of ways. For example, it 
produces a well-known inf-sup compatibility condition for mixed approximations 
that makes it difficult to achieve high accuracy with simple kinds of discretization. 
On the other hand, the incompressibility constraint is responsible for the energy 
inequality, an estimate which is fundamental to global existence theory. 

In this article we study global existence and stability questions for a non-de- 
generate parabolic system that is known to extend Navier-Stokes dynamics in a 
well-posed manner to velocity fields that are not necessarily divergence-free. This 
system appeared recently in [121 , and begins to explain the good performance of cer- 
tain numerical schemes where the pressure is computed by solving boundary-value 
problems [M] . The idea to determine pressure by solving boundary- value problems 
was also a feature of an earlier analytical study by Grubb and Solonnikov [6l[7] , and 
the system we consider is equivalent to one of their several 'reduced' models. 
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Explicitly, we study the initial-boundary value problem 

(1.1) 9tu + M • Vu + Vp - Au = in rj, 

(1.2) Vp= (/-P)(Au- VV-M-u- Vw), 

(1.3) u = on dn, 

(1.4) u ^ uq in f2, when t = 0. 

Here w = m(x, t) is the velocity field, p = p{x, t) the pressure, and P is the standard 
Leray projection of L^(r2,R'^) onto the subspace of divergence-free vector fields 
which are tangential at the boundary. For simplicity we have taken the kinematic 
viscosity to be unity and omitted body forces. 

For the system (|1.1|) - ()1.4|) . neither the initial data nor the solution are required 
to be divergence free. Equation (|1.2p defines the pressure gradient, and replaces 
the incompressibility constraint 

(1.5) V-M = inVL, 

that appears in the standard incompressible Navier-Stokes system. However, if 
initially V • uq = 0, then the incompressibility constraint (jl.Sp holds for all time. 
This follows because (ll.ip - dl.Sp show that V • u satisfies the heat equation with 
no- flux boundary conditions: 

{dtV ■ u = AV • u in f7, 
d 
— V-u = for a; e arj,i > 0, 

ov 

where ^ denotes the derivative with respect to the outward unit normal to 
Thus if V • Wo = 0, then the system (|l.ip -- (|1.4p reduces to the standard incompress- 
ible Navier-Stokes equations, and in this sense we say that the system p.ip - (ll.4p 
extends the dynamics of the standard incompressible Navier-Stokes equations. 

Of course, the dynamics of the standard incompressible Navier-Stokes equations 
could alternately be extended by completely omitting the VV • u term from (|1.2p . 
However, the presence of this term is crucial to the theory for two reasons. First, 
the AV • u term in (|1.6p is a direct result of the VV • u term in (|1.2p . and provides 
exponential stability of the divergence free subspace. Thus, from a numerical per- 
spective, errors in the divergence will be exponentially damped. The second, and 
perhaps deeper reason, is the essential role played by VV • u in the well-posedness 
results of [B1[71[T2]. We elaborate on this below. 

Define the Stokes pressure gradient, \Ipg{u), by 

(1.7) Vp,(u) = (/-P)(Am- VV-u). 

In context we often use pg to denote Ps{u)- Given any u £ H^{^, M''), the function 
Ps (u) is determined as the unique mean-zero solution to the boundary- value problem 

(1.8) Aps =0 inn, V -Vps ^ V ■ {A~ VV-)w on dVL. 

Without the VV • u term, the boundary condition in (jl.Sp would not make sense 
for all u G H^{n). With the VV • u term, however, Au — VV ■ u is and 
divergence- free; hence a standard trace theorem [U Proposition 1.4] makes sense 
of the boundary condition in H~^^'^{dfl). The Grubb-Solonnikov [6l[7] approach is 
based on using the boundary-value problem (|1.8p to determine the contribution of 



^Note that V ■ lio = Jqq uq ■ u = 0, and so the compatibihty condition for 1 11.611 is satisfied. 
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\7ps{u) to Vp, and proves well-posedness of (|l.ip - (|1.4|) using a theory of parabolic 
pseudo-differential initial-boundary value problems in L^'-based Sobolev spaces. 

The convergence arguments of Liu et al. |12) . on the other hand, result in a 
comparatively simple local well-posedness proof for (|l.ip - (|1.4p for initial velocity 
in HQ{n). This proof is based instead on the expression of the Stokes pressure 
gradient as a Laplace-Leray commutator: 

(1.9) Vps = {AP - PA)u. 
This follows directly from p. 71) using the fact that 

(1.10) W-u = A{I-P)u. 

Even in this approach, the 'extra' VV -w term in (|1.7p is directly responsible for the 
commutator representation (|1.9p . The key idea used in [H] is to treat the Stokes 
pressure gradient as the Laplace-Leray commutator (|1.9p . and show (Theorem l3.11 
below) that it is dominated by the Laplacian (cf. p.ip ) to leading order. 

While the methods of [6l[7l[T2] effectively address local well-posedness of (|l.ip - 
(|1.4p , they do not address global existence or stability. For the standard incompress- 
ible Navier-Stokes equations in three space dimensions, global existence of strong 
solutions is a well-known fundamental open problem [SlE]. However, classical re- 
sults establish global existence and regularity if if the flow is two-dimensional [TT] , 
or the initial data is suitably small [4l[TT]. 

In this paper, we establish a few such global existence results for the system (|l.ip - 
(|1.4|) . The main difficulty in proving a small-data global existence result for (|l.ip - 
(|1.4p is not the nonlinearity. The root of the problem is that the linear terms are 
not coercive under the standard inner product. We remedy this difficulty by 
using the commutator estimate in [T^ to construct an adjusted inner product under 
which the linear terms are coercive. This allows us to establish global existence for 
small initial data in two or three dimensions, and unconditional global stability 
of a time discrete scheme for the linear equations. This leads to an improved 
understanding of how the divergence and pressure can be handled to obtain energy 
estimates capable of controlling the nonlinear terms. 

In 2D, we can extend our small-data global existence results to initial data with 
small divergence. For arbitrary initial data, we can add a sufficiently large diver- 
gence damping term to (|l.ip - (|1.4p to obtain global existence. However, presently 
we are not able to prove global existence for (|l.ip - (|1.4p for arbitrary initial data. 
The difficulty is that for the energy balance using the standard inner-product, 
the non-linear term is skew-symmetric, and does not contribute; however, the linear 
terms are not coercive. On the other hand, for the energy balance using the adjusted 
inner products we consider, the linear terms are coercive; however, nonlinearity is 
no longer skew symmetric, and contributes non-trivially. 

Coercivity of the linear terms (albeit under a non-standard inner product) al- 
lows one to treat (|l.ip - (|1.4p as a non-degenerate parabolic system. While this has 
helped simplify existence theory and the analysis of certain numerical approxima- 
tion schemes, some other questions apparently become more difficult. In particular, 
while global existence of the standard incompressible Navier-Stokes equations is well 
known in 2D, it remains open for (|l.ip - p.4|) for general (2D) initial data. 
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2. Main results 

2.1. Coercivity of the extended Stokes operator. In the study of parabolic 
problems, an extremely useful (and often crucial) property is coercivity of the un- 
derlying linear operator. For (II. ip . the linear operator in question is the extended 
Stokes operator, A, defined by 

(2.1) Am = -Au + Vps(w) = -PA-it- VV-u. 

Note that the last equality follows from the identity (jl.lOp . Under periodic bound- 
ary conditions, the extended Stokes operator A is coercive. Indeed, under periodic 
boundary conditions, PA ~ AP, and so 

(2.2) {u,Au)^\\Vu\\l,, 

where (•, •) denotes the standard L'^ inner product on the torus. 

Under no-slip (0-Dirichlet) boundary conditions, the situation is surprisingly 
more complicated. The extended Stokes operator fails to be positive, let alone 
coercive, under the standard inner product. To briefiy explain why, observe 
that for ueH^n H^{n, R''), 

(2.3) {u,Au)= / u-Au^ / |V-it|^ + / u-Vps. 

Jn Jn Jn 

Now if V • u 7^ 0, the second term on the right need not vanish. In view of the 
commutator relation (|1.9|) . one might expect ||Vps|li2 to be dominated by ||Vm||^2- 
This, however, is known to be false, and control of the Stokes pressure Ps requires 
more than one derivative on u. Consequently, if V • m ^ 0, then the second term 
on the right of (|2.3p can dominate the first, and destroy positivity of A. 

Since our primary interest in the extended Stokes operator is to study (|l.ip - (|1.4l) . 
and the divergence of solutions to (|l.ip - (|1.4p is well controlled, one may hope to 
rectify non-positivity of A by a coercivity estimate of the form 

(2.4) {u,Au)^e\\Vu\\l,~C\\V -uWl,. 
But again, this turns out to be false. 

Proposition 2.1 (Failure of Coercivity). Let f2 C be a bounded, simply con- 
nected domain. For any £,C ^ 0, there exists a function u G C'^{Vl) such 
that 

(2.5) u = Qondn, and {u, Au) e\\Vu\\\., ~ C\\\7 ■ u\\\.,. 

The key idea in the proof is to identify the harmonic conjugate of the Stokes 
pressure as the harmonic extension of the vorticity. Since this is independent of 
our main focus, we present the proof of Proposition 12.11 in Appendix \^ towards 
the end of this paper. We remark, however, that if u G n H, where 

H = {v e L^iQ) \ v = Pv} = {v e L^{VL) I V • w = 0, and w • = on 9^}, 

then the second equality in (|2.1I) shows that the extended Stokes operator A reduces 
to the standard Stokes operator —PA. In the space n Hq n H coercivity of the 
standard Stokes operator is well known. Namely (|2.2I) holds for all u e iJ^niJg ni? 
(see for instance [H Chapter 4]). Unfortunately, when we consider vector fields for 
which u ^ H, Proposition 12.11 shows that coercivity fails for the extended Stokes 
operator. 
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The key to global existence results for the nonlinear system (|l.ip - ()1.4p is to 
remedy the negative results in Proposition ETT] in a manner that interacts well with 
the nonlinear term. This can be done by introducing a stabilizing higher order 
term, and a compensating gradient projection term, as we now describe. 

For any u G H^{fl) define Q{u), the primitive of the gradient projection, to be 
the unique mean zero function such that 

VQ{u) = {I- P)u. 

Given constants e, C > 0, we define an iJ^ -equivalent inner product ((•, •))^ ^ by 

(2.6) ((u, v))^ c = ^) + ^(Vu, Vv) + C{Q{u),Q{v)), 

where (•,•) denotes the standard inner product on L^{fl). Our main result shows 
that for all e sufficiently small, we can find C large enough to ensure coercivity of 
A under the inner product {{■, ■))^ q. 

Proposition 2.2 (i?^- equivalent coercivity). Let C K.'' be a domain. There 
exists positive constants Sq = So{fl) and c — c{fl) such that for any e G (0, eg), there 
exists a constant = C^{^1) > 0, such that for the inner product (■, ■)^ defined by 

we have 

(2.7) {u,Au)^ > i [\\Wu\\l, + e\\Au\\l, + C^VqWl.) 

for all u G H H^. Consequently, there exists a constant C'^ = C'^{e, fi) such that 

(2.8) {u,Au)^^ -{u,u)^, and (u, Au)^ ^ -i^(Vu, Vu)^, 

c 

for all ueH^nH^. 

We prove this Proposition in Section [3l The main ingredient in the proof is an 
estimate for the Laplace-Leray commutator (|1.9p that is proved in [12] and stated 
in Theorem 13 . 1 1 below. A couple of further consequences of this Theorem are worth 
mentioning here. First, A is invertible on with compact resolvent (Lemma l3.3l) . 
And, due to Theorem 13.11 and the self-adjointness of the Laplacian, an elementary 
result about sectorial operators [H Theorem 1.3.2] directly implies that A is a 
sectorial operator on with domain D{A) = £)(— A) = H Hq. 

The result of Proposition 12.21 raises the question of whether coercivity of A can 
be obtained in a space with less regularity than by using an equivalent inner 
product. In this regard we have two remarks. First, in Proposition 13.51 we will 
describe an inner product (•, •)^ for which A is coercive that is equivalent to the 
usual inner product on the space 

Hdw = {v e L'^{n) I V • t; G and v v = on dQ}. 

Second, we expect that a bilinear form defined by 

(2.9) {u,vy:'^{A-'/'u,A-'/\)^ 

determines an L^-equivalent inner product under which A is coercive. Coercivity 
for u G D{A) would follow from Proposition 12.21 and continuity by well-known 
interpolation estimates. However, an L^-coercivity bound {u,u)'^ > c|jM||^2 appears 
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not to be easy to prove — it may involve proving A has bounded imaginary powers 
(see [2]) in order to establish the expected characterization D{A^^^) = H^. 

In any case, unfortunately the inner products (■, ■}^ and (•, do not seem to in- 
teract well with the nonlinearity in (jl.ip . Thus for questions of global existence and 
stability for the nonlinear extended Navier-Stokes equations and their discretiza- 
tions, it is more convenient to use the inner product in Proposition 12.21 The rest 
of the paper can be read independently of Proposition 13 . 51 or its proof. 

2.2. Energy decay for the extended Stokes equations. A first step to global 
existence results for (ll.ll) - (|1.4l) . is the study of long time behaviour for the under- 
lying linear equations. These are the extended Stokes equations: 



(2.10) 



' dfU- Au + \7ps{u) ^0 inn, 

u{x, 0) = uo{x) for X e f2. 



A direct consequence of Proposition [2lT] is that the energy of solutions to (|2.10p can 
increase, at least initially. 

Corollary 2.3. There exists uq G C^(S7) with uo — on dfl, and tQ > such that 
the solution u to p.lOp with initial data uq satisfies 

\Hto)\\L'2 > \\u0h2. 

The proof of Corollary 12 .31 can be found at the end of Appendix \A\ following the 
proof of Proposition 12.21 

In contrast to the extended Stokes equations, solutions to the standard Stokes 
equations (with initial data in H) always have monotonically decaying norm. 
This follows because if u{t) € H, then multiplication by u and integration by parts 
produces the standard energy inequality 

(2.11) ld,\\uml. + \\Vu\\l.=0. 
The Poincare inequality now yields strict exponential decay 

(2.12) \\u{t)\\l,^e-^'\\uo\\l.. 

for all solutions to the standard Stokes equations with initial data in H. 

Despite the counter-intuitive initial energy increase, the extended Stokes system 
is a well-posed, non-degenerate parabolic system. This was proved in f6l!T2], and 
is a direct consequence of Theorem 13.11 Indeed, since A is sectorial it generates 
an analytic semigroup e^"^*, showing well-posedness of the initial-boundary- value 
problem (|2.10|) . Because no eigenvalue of A has non-positive real part by Propo- 
sition 12.21 one can quickly show that while the energy of solutions to (j2.10|) 
can increase initially, it must eventually decay exponentially. Explicitly, this means 
that solutions to (|2.10p must satisfy 

(2.13) \\u{t)\\l,^Ce-^'\\u4l, 

for some constants C, c > 0. 

To digress briefly, we remark that with a little work, one can explicitly charac- 
terize the spectrum of A. Indeed, if As denotes the (standard) Stokes operator 
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with no-slip boundary conditions, and An denotes the Laplace operator with ho- 
mogeneous Neumann boundary conditions, then 

a{A)^aiAs)Uc7{~AN)~{0}. 

Seeing cr(A) is contained in the right hand side above is immediate. The reverse 
inclusion requires a little work, and was communicated to us by Kelliher jl3| . 

Unfortunately, an abstract spectral-theoretic proof of p.lSp is not of direct help 
for studying the stability of time-discrete schemes, which was a primary motivation 
for introducing these equations. Further, p.lSp does not recover ()2.11|) for solutions 
with initial data in H. For this reason, we search for a direct energy-method proof 
of (|2.13p . and for an idea which also allows the study of time discrete schemes. 

Observe first that if we multiply (I2.10p by u, integrate, use the commutator 
estimate (13. ip and Gronwall's lemma, we obtain exponential growth, not decay, of 
||u||^2. If we involve a higher derivative, coercivity of A in Proposition 12.21 (or 
Proposition 13. 5|) and Gronwall's lemma guarantee eventual exponential decay of 
ll'^llifi (o'" ll^llffd- Ho'wfiver, for (j2.10|) . we can obtain a more satisfactory decay 
estimate by considering non-quadratic form energies. 

Proposition 2.4. Let u be a solution to (|2.10p with uq £ H^{il). Then for any 
e > 0, there exists constants ci — ci(il) and C2 — C2(ri,e), such that ci,C2 > and 

(2.14) dt£c^^a,{u) + SU^) ^ 
where £ci,c2 '"^'^ ^'e o'^g defined by 

(2.15) f,„,,(u) IM\\u\\l.+c,\\Vu\\^4\/Q{u)\\^.+C2\\VQ{u)\\l., 

(2.16) £Uu) '^{2 e)\\Vu\\l, + \\Au\\^4VQiu)\\^, + \\AQ{u)\\l,. 

The proof of Proposition[2]4]is in Section|4l While (I2.14p docs not imply eventual 
exponential decay controlled only by the norm as in ()2.13p . it does provide an 
estimate that reduces to the energy inequality for extended Stokes equations (|2.1ip 
when the initial data is in H. To see this, note that if uq E H, then (/ — P)u{t) — 
for alH > because V ■ u satisfies the heat equation (11.61) . Consequently Q{u) = 0, 
and equation (|2.14p reduces to 

dt\\u\\l, + {2-e)\\\/u\\l,^Q. 

Thus in the limit e — ?> 0, we naturally recover the energy decay for the Stokes 
equation (equation (|2.1ip ) for initial data in H. 

We also notice that the 'energy' £ci,c2 of solutions must in fact decrease expo- 
nentially. This is because /q Vw = = J^Q{u), and so the Poincare inequality 
can be apphed to both the terms ||Vu||^2 and ||Q(u)||^2. Thus equation ()2.14p 
immediately implies 

for some small constant c = c{ci,C2,e,Q). Unfortunately, however, for the ex- 
tended Navier-Stokes equations, the 'energy' £ci,c2 does not interact well with the 
nonlinearity. 

2.3. Uniform stability for a time-discrete scheme. Before moving on to the 
non-linear system (|l.ip - (|1.4p . we study stability of a time-discrete scheme for (j2.10p . 
of the type treated in [12 . One main motivation for studying the system l|l.ip - (|1.4p . 
or the linear system (|2.10p . is that this kind of time-discrete scheme is naturally 
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implicit only in the viscosity term, and explicit in the pressure. We will show 
that the ideas used in the proof of Proposition 12.21 give globally uniform stability 
estimates for such time-discrete schemes. 

Given an approximation u" to the velocity at time n 6t, we determine Vp" from 
the weak-form Poisson equation 

(2.17) (Vp",V(^) = (Au"-VV •-«" + /", V(^) VipeH\n). 
Now, we determine by solving the elliptic boundary value problem 

( - It" 

(2.18) <^ St 
[ = on on. 

where /" — ^ Inst'^'^^* ^® ^ time-discretized forcing term. 

Proposition 2.5. Let fl be a bounded domain in M'^, d — 2,3, with boundary. 
Then there exist positive constants Kq, e^Ce^CyC , depending only on il, such that 
whenever < dt < kq, then for all N > we have 

(2.19) \\u^\\l.+e\\S7u^\\l+C,\\Vq''\\l.+ 



Au"+i + Vp" = /" in n, 



1 ^ 



C 

k=0 



lL2 



St 



^\\u%.+e\\VuX.+C,\\\7q''\\l. 



N 



+ CSt[\\Vu"\\l.+e\\Au°\\l.+Y,\\f'^\\h 



k=0 



an 



d 



(2.20) h^||i.+e||V^.^||i.+Ce||Vg^||i. 



^ (1 - CStr (\\uX. + eW^uX. + a||VgO||i. 



+ c'J2\\fk\\U^-cstf-'-''st. 

k=0 

The proof of this proposition is in Section [5l 

2.4. Global existence results for the extended Navier-Stokes equations. 

When one seeks an energy estimate for (|l.ip - p.4p . multiplying (|l.ip by u, the 
nonlinearity produces the term 

(2.21) / u • (u • V)u = /" (V • ^i)|w|^ 

In general this is non-zero, but is morally harmless since V ■ m is a solution of 
(|1.6p and is well controlled. This is indeed the case in two dimensions, but under 
periodic boundary conditions (see Proposition 12.91 and the remark following it). 
The key ingredient for proving global existence for periodic boundary conditions 
is the coercivity (|2.2I) of the linear terms. Consequently, despite the extra non- 
linear term arising from (|2.2ip . the energy balance closes and the well-known 
existence results for the standard incompressible Navier-Stokes equations continue 
to hold with minor modifications. 



COERCIVITY AND STABILITY FOR EXTENDED NAVIER-STOKES 



9 



The situation is more complicated mider no-slip bomidary conditions, however, 
since now coercivity (|2.2|) fails. To get any mileage from the linear terms, we need 
to use an inner-product under which the linear terms are coercive. Using the inner 
product in Proposition 12. 2[ and a 'brutal' estimate on the nonlinearity, we can 
obtain a two or three dimensional small-data global existence result. 

Theorem 2.6 (Small data global existence). Let d = 2 or 3, il C M'^ be a bounded 
domain with boundary. There exists a small constant Vq — Vo{fl) > such that 
if uq £ HQ{n) with 

II^oIIhi < Vo 

then there exists a global strong solution to - (11.41) with 

(2.22) u e L^{0,T; H^{n) n H^{n)) n H\0,T; L'^{n)). 

for any T > 0. Consequently u G C([0, oo); Hq) and V ■ u e C°°{{0, oo) xfl). 

The proof of this theorem is in Section [6l Two-dimensional global existence, 
however, poses a different problem. A key ingredient in 2D global existence for 
the standard incompressible Navier-Stokes equations is the energy balance: the 
nonlinearity cancels, and doesn't contribute! Unfortunately, for (|l.ip - (|1.4p . the 
L^-energy balance doesn't close because of the higher order contribution from the 
Stokes pressure gradient. 

In the absence of an energy inequality, we are only able to prove a pertur- 
bative result. If the initial data is divergence free, then (ll.ip - (|1.4p reduces to the 
standard incompressible Navier-Stokes equations, for which 2D global existence is 
well known. Thus for initial data with small divergence, we can prove 2D global 
existence for (fLTjl -pn i) . 

Theorem 2.7 (Small divergence global existence in 2D). Let C be a bounded 
domain, vq G i^o(f^) with V ■ vq — be arbitrary. There exists a small constant 
Uq — Uoifl, \\vq\\ ffUQ-f) > such that if 

(2.23) uoeHgin), PoUo = vo and || V • Mo||^2(f2) < J7o 

then there exists a global strong solution to (|l.ip - p.4p with initial data uq such 
that ^jm^i holds for allT >{). 

The operator Pq above is the i/Q-orthogonal projection of Hq{VI) onto the sub- 
space of divergence free vector fields, and is described in Section [7] along with the 
proof of Theorem l2.7l One strategy to avoid the small divergence assumption is to 
further damp the divergence. Namely, for arbitrary initial data (in 2D), if we add 
a strong enough divergence-damping term to (|l.ip - (|1.2p . we can guarantee global 
existence. 

Corollary 2.8 (Divergence-damped global existence in 2D). Let <Z M? he a 
, bounded domain and uq G H^ifl) be arbitrary. There exists a constant ag = 
C(o{flj II V • uo|lL2(f2')) > such that if a ^ uq then the system 

' dtu + P{{u ■ V)u) + Au + a{L - P)u = m n 

(2.24) \ u(x,t)=0 forxedn,t>0 

u{x,0) = uo{x), 

has a global strong solution u such that p.22p holds for all T > 0. 
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The main idea in proving Corollarv l2.8l is to verify that the divergence-damped 
extended Stokes operator Ba defined by 

(2.25) = A + a{I - P) 

is coercive, with coercivity constant independent of a. Consequently, the proofs of 
Theorems 12 .61 and 12 . 71 work verbatim for the system (|2.24[) . with constants indepen- 
dent of a. Combining these existence theorems, and using the added divergence 
damping gives Corohary 12. 8[ a better existence result as an easy corollary. We 
devote Section [8] to the coercivity of Ba (Proposition [831, ^'^'^ the proof of Corol- 
lary El 

So far, our two-dimensional global existence results under no-slip boundary con- 
ditions required either a small initial divergence assumption, or an additional strong 
divergence damping term. Such requirements are not needed under periodic bound- 
ary conditions, primarily because of (|2.2p . We observe, then, that the identity (|2.2p 
will still hold in domains with boundary, provided we consider functions u with 
boundary conditions 

(2.26) Pu • r = on and m ■ = on dVL, 

where v and r arc the unit normal and tangential vectors respectively. These 
boundary conditions (|2.26p reduce to the usual no-slip conditions in the physically 
relevant situation where u = Pu. 

Armed with (12. 2p . we obtain a 2D global existence result without a smallness 
assumption, or any additional divergence damping. 

Proposition 2.9. Let fl C M.^ be locally Lipschitz and bounded, and let uq G 
H^{fl). There exists a time-global strong solution to (|l.ip - (|1.2p with initial data 
Uq and boundary conditions p.26p . 

We prove the identity l|2.2p and Proposition 12.91 in Section [HI The proof of 
Proposition 12.91 emphasizes another (analytical) advantage of the boundary con- 
ditions (|2.26p . Under all the boundary conditions we consider (no-slip, periodic, 
and (j2.26l) ) the evolution equation for the gradient projection is always linear, self 
contained, and decays at an explicitly known rate. The evolution equation for the 
Leray projection (equation (19. 2p ). is coupled to the gradient projection; however 
the coupling terms are harmless. What causes trouble under the no-slip bound- 
ary conditions is that the evolution of the Leray projection is also coupled to the 
gradient projection through boundary conditions! This proves problematic in the 
case of 2D global existence. On the other hand, periodic boundary conditions, or 
the boundary conditions (|2.26p provide an explicit de- coupled boundary condition 
for the Leray projection, which simplifies the analysis greatly. Unfortunately, the 
price paid is that the boundary conditions (|2.26p are much harder to implement 
numerically. 

3. Coercivity of the extended Stokes operator. 

As mentioned earlier, the extended Stokes operator is not coercive under the 
standard inner product. However, it is coercive under a non-standard, but H^- 
equivalent, inner product. This is the main tool we use in studying the extended 
Navier-Stokes. The aim of this section is to prove Proposition l2.2l (coercivity under 
the adjusted inner product). The main ingredient in the proof is the following 
estimate on the Laplace-Leray commutator. 
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Theorem 3.1 (Liu, Liu, Pego [12 ). Let be a connected, bounded domain with 
boundary. For any S > there exists Cs ^ such that 

(3.1) \\yps{u)\\i, ^il^^) ii^^iii^ + ^^ii^^iii^ 

for all u e H^nH^ifl). 

We refer the reader to [T^] for the proof of Theorem 13.11 

3.1. _ff ^equivalent coercivity on D{A'^). The idea behind the proof of Proposi- 
tion 12.21 is to use Theorem 13.11 and prove coercivity assuming the 'extra' boundary 
condition An G Hq. We will later use an approximation argument to prove the 
Proposition for all fl Hq functions. 

Lemma 3.2. For any e > sufficiently small, there exists a constant c — c{n), 
independent of e, and a constant C'^ = Ce(r2) > 0, depending on e and fl, such 
that ((2J)) holds for all u e H'^ n Hq, such that Au e Hq. 

Proof. Observe first that there exists a constant C = C{il) such that for all u G 
iJ^ n H^ we have 

(3.2) \\Vps{u)\\^.^C\\Au\\^.. 

While this immediately follows from Theorem 13.11 and the Poincare inequality, we 
can see it directly from (II. 7|) because, due to elliptic regularity, 

(3.3) \\^Ps{u)\\l. ^ C (llA^lli. + II VV • uWl^) ^ C||A«||i.. 

Now let u e H"^ n Hi be such that An e H^, and q = Q{u) be the unique mean 
zero function such that Vq — {I — P)u. Then 

(3.4) (m, Au) = {u, -Am) + (u, Vp,) = || Vu||^2 + (Vg, Vp,) 

> llV^lli, -||Vg||^.||Vp,||^. ^ llVulli. -C||Vq||^.||A7.||^. 

^ \\Vu\\l.-^^\\Au\\l,-C,\\Vq\\l. 

where the second last inequality followed from (|3.2I) . and is some constant de- 
pending only on 17 and e. 

Since Au = on dVl by assumption, we can integrate the _ff^-term by parts. 
This gives 

(3.5) (Vu,V^u) ^ ~{Au,-Au + Vps) > ^||Au||^2 - Ci||Vuf 

8 

where Ci is the constant that arises from Theorem 13. II 
Thus if £ < 2^ ; equations p.4p and p.Sp give 

(3.6) {u,Au)+e{Wu,WAu) > ^\\\/u\\l, + ^\\Au\\l, - C,\\Vq\\l. 

Now let r be the unique mean zero function such that Vr = (/ — P)Au. Observe 
that 

(/ - P)Au = (/ - P){-PAu - VV ■ u) = -VV • u. 
Since V • m = 0, we must have r = —V • u = —Aq. Thus 

{q,r)^{q,-Aq) = \\Vq\\l,. 
Combining this with (|3.6p we get (|2.7p as desired. □ 
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3.2. Properties of the extended Stokes operator. Consider the extended 
Stokes operator A as an operator from into L'^{n) with domain D{A) = 

n iJp. In this context, we recall that Proposition 12.21 asserts (|2.7p for all 
u G D{A)\ however, Lemma [3.21 only proves (j2.7p for all u e D{A^). To address 
this gap, and finish the proof of Proposition 12.21 we need a few basic properties of 
the extended Stokes operator. 

Lemma 3.3 (Regularity and invertibility) . The extended Stokes operator A has a 
compact inverse. Furthermore, there exists a constant c — c(f2) > such that 

(3.7) ^I|u|Ih2 =^ \\Au\\^^ s$ c||u||^2, for all u e Ci 

Proof. Our first step is to obtain estimates for the operator A + \I with A large 
enough. For an arbitrary u e H'^ n Hi, let f = {A + \I)u. Multiplying by -Au 
and integrating gives 

A||Vm||^2 + ||Au||^2 = / Vp,-Audx- / fAudx 
Jn Jn 

^ l\\Au\\l. + i||Vp,||i. + ^\\Au\\l. + 4||/||i. 

G + ? + i^) iiA"iii^ + ciiv«iii=+4|i/iii. 

where the last inequality followed from Theorem 13. ![ and c = c{fl) is a constant. 
This gives 

i\~c)\\yu\\l. + ^\\Au\\l.^i\\f\\l. 
Thus, when A > c, we immediately see 

(3.8) \\u\\^2^C\\f\\^2=C\\{A + XI)u\\^2, whenueH^nHl. 

One can use the last relation to check that A is closed. We claim further that 
A + XI is surjective for some large enough A. This can be proved by a Neumann- 
series perturbation argument based on the identity 

(3.9) A + XI = {I + B){XI ~ A), 
where B = Wps o (A/ - A)"^ That is, 

Bu = Vps{v), v = {XI-A)-^u. 

It suffices to prove that the operator norm of B on is strictly less than one, if A is 
positive and large enough. By easy energy estimates, we have that A||w||^2 < ||M||i2 
and ||Ati||^2 < I|w|Il2. Then due to Theorem 13. II and interpolation, we have 

115^11^2= / \Vps{v)\-' ^p\\Av\\l.+Cp\\v\\%<{p + CpX-^)\\u\\%, 
Jn 

and the coefficient on the right is less than 1 for A large enough. Thus / + i? is an 
isomorphism on L^, hence A + XI is surjective. 

Further, the Rellich-Kondrachov compact embedding theorem and the bound 
p.Sp imply that A+XI has compact inverse. Since we have shown that the resolvent 
of A contains at least one element with a compact inverse, the spectrum of A consists 
only of (isolated) eigenvalues, of finite multiplicity (see for instance O Theorem 
III. 6. 29]). Thus to prove invertibility of A, it suffices to show that is not an 
eigenvalue of A. 
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To see this, suppose u G D{A) is such that Au ~ 0. Then —PAu = VV • u e 
i^(r2). Since the range of the Leray projection (by definition) is orthogonal to 
gradients, we must have VV • u = PAu = 0, and hence V • u must be constant. 
Since J^V - u = Jg^ u-iy = 0, this forces V • u = 0. Thus u = Pu, and is orthogonal 
to gradients. Since Au — 0, we have 

= / u- Audx^ - / u- Audx+ / Pu ■ Vps{u)dx = ||Vw||^2 + 0, 
Jn Jo. Jn 

forcing u = 0. Hence is not an eigenvalue of A, and we conclude that A is 
invertible. 

It remains to establish (13.71) . The upper bound follows immediately from (j2.ip 
and (|3.2p . To prove the lower bound, observe first that boundedness of A^^ implies 



(3.10) \\A-'u\\^. ^ C\\u\\ 



L2J 



for some constant C — C{il), which we subsequently allow to change from line to 
line. Thus using the operator identity 

= (A + A/)-i (/ + AA-1) , 

and the inequalities (13.81) . p.lOp we see 

\\A-^u\\^, = \\{A + \I)-^{l + XA-^)u\\^, 

^C\\{l + XA-')u\\^,^C{l + CX)\\uh. 

proving the lower bound in p.7p . □ 



Lemma 3.4. For the extended Stokes operator, D{A^) is dense in D{A). 

Proof. Let u G D{A), and v = Au. Since v e i^(fi), we can find v„ G fl Hq 
such that (vn) v in L^. Since D{A) = n Hq by letting u„ = A^^Vn, we have 
Un G D{A'^). Finally, by Lemma [23] we see 

\\Un - u\\^2 < c\\AUn - Au\\ ^2 = c\\Vn - v\\^2 ^ 0, 

concluding the proof. □ 



3.3. iJ^-equivalent coercivity on D{A). Lemmas 13.21 and 13.41 quickly imply 
Proposition [221 



Proof of Proposition \2.2[ Let u G D{A). By Lemma [3.41 there exists a sequence 
Un G D{A'^) such that (u„) — > m in H^. By Lemma [3.2[ there exists constants 
c{n), Ceifl) > such that 

{Un,AUn), ^ i (||VM„|j^2 +e||Au„||^2 +Ce||V(7„||'2) , 

where qn is the unique, mean-zero function such that V(7n — (1 — P)un. Since 
(m„) u in H^, taking limits as n ^> oo yields (|2.7p . Now using the Poincare 
inequality, (|2.8I) follows. □ 
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3.4. Hdiv equivalent coercivity. We conclude this section by proving coercivity 
under an i/div-cquivalcnt inner product. The rest of this paper is independent of 
this result and its proof. 

Proposition 3.5 (i/div-coercivity). Let C be a domain. There exists 
positive constants Sq > and c = c(f2) such that for all e £ (0,eo), there exists a 
constant — Cf^{Q) > such that the following hold. 

(1) Let (•, •)^ be defined by 

(3.11) (u, v)'^ = {u, v) + £(V • u, V • w) + Ce{A-\, A-\)^ ~ 

-{u,VVs{A-\))-{^Vs{A~^u),v), 

where {u,v) is the standard inner product on L'^{fl), and i-,-)^ denotes the 
inner product from Proposition \2.Si Then 

(3.12) -Juh. + £|| V • u\\l, < {u, u)[ ^ c(l + CMAl. + £|| V • u\\l. 

for any u £ L^(i7) with V • u £ i^(r2) and u ■ v — on dfl. 

(2) For any u £ H with V • u £ , we have 

(3.13) ||Vu||^. + |||VV.u||^. 

(3.14) and {u,Au)' > —{u,u)' . 

c 

Proof. We begin by proving p.l2p for all u £ D{A). Density of D{A) in and 
a standard approximation argument will now establish p.l2p for all u £ -ffdiv We 
will assume and c are constants that can change from line to line, provided their 
dependence on parameters is as required in the Proposition. Let u £ D{A), and 
V = A~^u. Then from (|3.1ip we have 

{u, u% - {Av, Av) - 2{Av, Vpsiv)) + e||V • u\\l, + C,{v, v)^ 

= {-Av + Vp,{v), -Av - Vpsiv)) + e\\V ■ u\\l2 + Ce{v,v), 

^^■^^^ - II AHli. - \\yps{v)\\l. + e||V • u\\l. + CAv, v)^ 

^ ^\\Av\\%-c\\v\\%+s\\V-u\\l,+C,{v,v)^ 

where the last inequality followed from Theorem 13.11 and interpolation. Now since 
Av = u, we immediately see 

\\u\\l,_ = \\-PAv-^V-v\\l,^c\\Av\\l,. 

Finally, by definition of {v,v)^, we have {v,v)^ ^ II'^^IIl^- Thus if > c, the lower 
bound in equation (|3.12p will hold for all e > 0. 

For the upper bound in (|3.12p . observe that by definition of Vps, and Lemma l3.3l 
we have 

WVpsA-^uW^, s$ c\\A-^u\\jj-2 sc c||m||^2. 

Combined with the estimate ||^"^u||^2 ^ c||u||^2, which is also a consequence of 
Lemma 13.31 we immediately obtain the upper bound in p.l2p . 

Finally, it remains to prove the inequality p.l3p . We will prove (I3.13P for u £ 
D{A'^); since Z)(A^) is dense in D{A), the same approximation argument from the 
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proof of Proposition 12.21 will show that p.l3p holds on D{A). In keeping with the 
above notation, we again set v — A^^u. This gives 

(u, A-it)^ = (u.Au) - {u, Vps{u)) - {Vps{v),Au) + e(V ■ u,V ■ Au) + Ce{v,Av)^. 

We deal with the terms on the right individually. Combining the first two terms, 
the dangerous term involving the Stokes pressure cancels. This gives 

{u,Au) - {u,Vps{u)) = {u,~Au) = ||Vu||^2. 

For the third term, 

^{Vp,{v),Au) = -{^ps{v),{l - P)Au) - (Vps(f),VV-u). 

For the fourth term, observe that if u e D{A^), then Au = on dfl. Thus 
integrating by parts gives 

{y -u^y -Au) = -{vy ■u,Au) = ||vv-u||^2. 

Combining these identities we have 

(u, Au)[ = WVuWl, + {Vpsiv), VV • w) + e||VV • u\\l, + C,{v, Av)^ 

> llVwIli. + |||VV • u\\l, - l-\\yp,{v)\\l. + C,{v, Av)^ 

> llVz^lli. + |||VV • u\\l. - -^\\I^v\\l. + C,{v, Av)^ 

> llVz^lli. + |||VV • + - ||At;||i. + ^i^v^Av)^ 

where the last inequality followed from Proposition 12.21 and c\ is the constant 
in (12. 7p . We note that Proposition 12 . 21 guarantees that c\ is independent of £. Now 
we choose large enough so that — | > 1, giving 

(3.16) («, ^ llV^lli. + |||VV • uf^. + ||Az;||i. + ^(z;, Az;),, 

from which inequality p.lSp follows. 

Finally for p. 141) . observe that from p.lSp we have 

^C2||A«||^. +e||V-u||^. +Ce(f,v)„ 

for some constant C2 = ci(^Y Now using Proposition 12.21 and p.l6p . the inequal- 
ity ([3ll| follows. □ 

4. Decay of a non-quadratic form energy. 

This section comprises the proof of Proposition 12.41 addressing the long time 
behaviour of solutions to the extended Stokes equations p.lOp . The result and 
proof are independent of the rest of this paper. 

Proof of Proposition \2.4\ In this proof, we use C to denote an intermediate constant 
that depends only on Q. whose value can change from line to line. We use Ci, C2, . . . 
to denote fixed positive constants that depend only on Jl, whose values do not 
change from line to line. 

As usual, let q = Q{u) be the unique mean zero function such that Vg = {I—P)u. 
We begin by establishing the energy inequalities 

(4.1) d4u\\l. + 2\\Vu\\l. ^C4yq\\^.\\Au\\^., 
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(4.2) dt\\Vu\\l.^\\\^u\\l.^C\\Vu\\l,, 

(4.3) dt\\Vq\\l,+2\\V -uWl.^Q. 

Before proving the above inequalities, we remark that the balance (j4.ip does 
not close by itself. On the other hand, the balance (|4.2p closes, but does not 
give decay. A combination of the norms, however, gives us the desired exponential 
decay. 

For the proof of (|4.ip . multiply (|2.10p by u and integrate over 51 to obtain 

\dt\\u\\l, + \\Vu\\l, = - f u-Vps ^ - f Vq-Vps ^ Ci\\Vq\\^4Au\\^,, 
^ Jn Jn 

where we used (|3.3p . This establishes (14. ip . 

Turning to (|4.2p . we multiply ()2.10p by —Au and integrate over to obtain 

(4.4) ldt\\Vu\\l, + \\Au\\l, = f Vps- Au. 

Using Theorem l3.11 we know that for any S > there exists a constant Cs — Cs{i^) 
such that 



Hence 



ypAh^ [^]\\Au\\l,+Cs\\Vu\\l,. 



(4.5) 



Vps • Au 



Choosing S — ^, equation (j4.4p reduces to (|4.2p as desired. 
Finally for gS]), we apply (/ - P) to (I^JU)) to get 

dtVq + (/ - P){-Au + Vpsiu)) = 0. 

Since V • m = Aq and 

(/ - P){-Au + Vpsiu)) = {I - P){PAu - VV • u) = VV ■ u = VAg, 

we see 

dtVq - AVq = 0, 

and hence dtq — Aq ~ C{t), where C is constant in space. Now, since u = on dil, 
we must have ^ = on 957. This means Aq — 0; since J^q — by our choice 
of q, we must have C(t) = 0. Thus we obtain 

d 

(4.6) dtq - Aq = in Q, with —q = for a; £ dQ. 

ov 

Multiplying by —Aq and integrating over 57 gives (|4.3p as desired. 

Now we combine ()4.ip - ()4.3p to obtain the desired exponential decay. First 
from (j4.2p . (|4.3p and the Poincare inequality we have 

5*l|Vg||i2 + -^||V(7||^. ^0, 

Thus 

at(||Vw|l^.I|VgI|^.) + -^||A^.|l^.||Vg|l^. ^C4l|V7.|l^.I|VgI|^.. 



COERCIVITY AND STABILITY FOR EXTENDED NAVIER-STOKES 



17 



where C4 = C3 — Using this in (|4.ip we see 

dt[\\u\\l,+C2{C, + l)\\Vu\\^4\7q\\^,) + 

+ 2||Vu||i. + ||A^.|l^.|lVg|l^. ^C5||V^.|l^.||Vq||i. 

where C5 = C2(Ci + 1)C4. Now letting ci — C2(Ci + 1), and C2 — 02(^2, e) to be 
chosen later, we see that 

dt (ll^lli^ +ci||Vu||^.||Vg||^. +C2||V<z||^.) + 

+ 2\\Vu\\l. + \\Au\\^. ||Vg||^. + 2c2|!V • u\\l. ^ e\\\7u\\l. + ' ""'^ 

where Ai is the best constant in the Poincare inequality 

X4Vq\\l.^\\Aq\\l.^\\V-u\\%. 
Thus choosing C2 = + i, we obtain (|2.14l) . □ 

5. Global stability of time discretization for the extended Stokes 

equations 

We devote this section to proving Proposition l2.5l The main idea again is similar: 
to introduce a stabilizing, high order term in the definition of the energies. 

Proof of Proposition \2.5i In the following we use Ci , C2 , . . . to denote fixed positive 
constants that depend only on CI, whose values do not change from line to line and a 
generic constant C whose value might change from one line to the next, depending 
only on f2. 

Let p" be the Stokes pressure for u" hence Vp" = (AP — PA)u" and thus 

[ Vp^' -Vip^ [ (Au" - VV • u") • V(^, Vip G H\n). 
Jo. Jo. 

Using (|2.17p with = and combining it with the last relation we obtain 

(5.1) w\\L^^w:\\L^ + \\n\L^- 

We derive first the discrete 77^ estimate just as in 12 . Taking the inner 
product of (Pl^ with -Au"+i gives 

^ (l|Vz."+ii. - ||V^"||i. + ||V^"+i - V^"||i.) + ||A^"+i||i. 
sc: ||Au"+i||^. (2||riL. + ||Vp«||^0 

^ |iiAu"+ii. + -iiriii. + i(iiA«"+ii. + iiv/;iii.) 

for all El > 0. This implies 

(5.2) 1 (||Vz."+ii. - l|V«"||i.) + ^l|V^"+i - Vh'I^. + 

+ (1 - ei)||Au«+i||i. -linii. + ||Vp';||i.. 
Fix any f3 G (5, §)• By Theorem 13. II we have 

||Vp:||i. <^/3||A^."||i. + ^||V^."|li.. 
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Using this in ()5.2|) and dividing by 2Cp we get 



(5-3) ( l|V^."+ii. - ||Vz."||i. + ||V^"+i - Vz.""^ 



2Cp5t V" ^ " llL2 - II V a 11^2 -r II V U - Vti 11^2 

<^iiriii^ + ^iiv«"iii., 

and we may assume that ei > is smaU enough so that 1 — £i — |/3>0. 

We continue by obtaining the discrete estimate. We dot the equation (|2.18p 
by m""*"^ in and obtain 
1 

+ ^l|V«"+^|li2 + I (||Vu"+l||i2 - ||V^"||i2) + \\\yU-\\l. 

if" - Vp") • - / (PD ■ u"+' - [ \/ps{u") ■ Vq 

< llrlL2||^^"+^IL2 + ||Vp,K)||^2||Vg"+i|1^2 

(5.4) ^ ^\\u"+X. + ^—^l-^WAu-Wl. + CiWnl + l|Vg"+ii2), 

where — and Aq is the principal eigenvalue of the Laplacian on 

with zero Dirichlet boundary conditions. 

Since p" = p'^ + by applying / - P to ([2381) we find that g" satisfies the 

time-discrete inhomogeneous heat equation 

(5.5) 1 ~ g") - Ag"+i = p^+' - p^. 

Then we find after testing with — Ag"+^ that, as above (and as in [121 page 1477]), 

(5.6) ^ fl|Vg"+i||i2 - ||Vg"||i2) + \\Aq"+'f ^ Ib^i -p';||i2. 



n+l 



St 

Also, since satisfies a Neumann boundary value problem, we have the estimate 

(5.7) ib^i -p^:\\i.. ^ ch"+i -«"iii2ii«"+^ -^"iii2. 

Now choose C2 large enough to ensure 

Ci||Vg"+i||i2«:^||Ag"+i||^2, 
and £2 small enough so that 4e2 < 1 ^ £1 — Combining (15. 6p and (15. 7p . we obtain 

(5.8) i(||Vg"+i||i2-||Vg"||i2)+||Ag"+i|r 



^^^l|A^"+i-A^"||i2+C3h'^ 



u 



^ -£^||A^"+i - A^"||i2 + C4||Vu"+i - V«"||i2 

for large enough constants C3 and C4. 

Assume that 5t is small enough so that {2Cf)C2Ci)5t < 1. Multiplying (|5.8p by 
C2, and adding it to (|5.3p and (|5.4p then gives 
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(5.9) ^ (||V<z"+ii. - l|Vg"||i.) + ^||Ag"+i||i. + ^ - \W1 

(liv."-||i.-||v."||i.) + i||v."||i. 



^2 /ii A ,n+l||2 11 A „n||2 \ , ( ^ , ^ \ n m||2 



Now summing from n = to A'^ the last inequality gives (for small enough St, and 
for a suitable constant C > 0) the claimed inequality p.l9p . 
We rearrange (|5.9p and obtain 

(5.10) j^(||Vg"+i||i.-(l-C<5t)||Vg"||i.) + 

- (1 - C6t)\\um.) + 1 (||VH"+ii. - (1 - C<50||V^"||i.) + 

+ (||Vu"+i||i.-(l-CTi)||V«"||i.) + 
+ (1 + 6t) (||Au"+i|li. - (1 - CTi)l|A^"|lL) + 
(||Ag"+ii. - (1 - (7<5t)||Ag"f ) +(1 - CSt)\\Aq-\\h ^ C\\r\\h 

provided that St is small enough, for suitable constants C and C. Defining 

a« = j^||Vg"||i. + l||u"||i. + l||V«"||i. 

+ ||Vu"||i. + (l + <5t)||Au"||i. + ||Ag"||i., 

()5.10p becomes 

(5.11) a„+l - (1 - CT^fln ^ CWfnWh. 

Solving this recurrence relation yields (|2.20p . □ 



6. Small data global existence for the extended Navier-Stokes 

equations 

This section is devoted to the proof of a long time, small data existence result 
(Theorem 12. 6p for the system (|l.ip - (|1.4p . As bounds for the linear terms have al- 
ready been established (Proposition 12. 2p . we begin with a bound on the nonlinear 
term. When obtaining energy estimates for solutions to ()l.ip . the explicit, expo- 
nential decay of V • it allows sharper estimates for many terms. However, in order 
to exploit coercivity of the linear terms, we are forced to use an i/^ -equivalent inner 
product. In this case, the 'worst' term that arises from the nonlinearity isn't aided 
by decay of V • u, and must be estimated brutally. Consequently, estimating the 
remaining terms similarly doesn't weaken the final result. Thus, we begin with a 
lemma that provides a 'brutal' estimate on the nonlinearity. 



20 



GAUTAM IYER, ROBERT L. PEGO, AND ARGHIR ZARNESCU 



Lemma 6.1. Let f,g,h E r\ H}^{Vl) with Q. C W^,d = 2,3 a hounded domain 
with boundary. Then there exists a constant C — C(f2) > such tha^ 

KP((/.V)g),/i)J =C|l/|l^,|lVg||^i||A/i||^.. 

Proof. Observe first that 

(P((/ ■ V)5), h)^ = (P((/ • V)g),h) + e(VP((/ • V)g), V/i), 

since (/ — P)P{{f ■ \7)g) ~ 0. Tlius to prove the lemma, it suffices to show the 
estimates 



(6.1) KW-V)5),/i)KC||/||^,||Vg||^i||/i 



lL2 



(6.2) and |(VP((/ • V)g), V/i)K C||/||^, ||V.g||^i ||A/i||^.. 

for some constant C — C{fl). 

The inequality (j6.ip follows directly from the Sobolev embedding theorem. In- 
deed, for any three functions /i, /2, /a, we kno^wQ 



(6.3) 



/1/2/3 



^C^II/l|lif=JI/2||^=.ll/3||^= 



provided ^ Si ^ 3, si + S2 + S3 ^ f and at least two of si, . . . , S3 are non-zero 
(see for instance the proof of Proposition 6.1 in [J]). Choosing si = 1, S2 = 1/2 
and S3 = 0, we have 

KP((/ • V)g),h)\ = |((/ . V)g,Ph)\ ^ CUWh. llVsll^i 

proving (|6.ip . 

For (j6.2p . we first integrate by parts and observe AP is a regular differential 
operator (identity (|1.10p ). Now we can integrate by parts again to obtain the 
desired estimate. Explicitly, 

(VP((/ • V)g), \7h) = -(AP((/ . V)g),h) 

- -((A - VV-) ((/ • V)g), h) = -((/ • V)5, (A - VV-) h), 

where all boundary integrals vanish because f,g,hE Hq. Now using (j6.3p with si — 
1, S2 — 1/2, S3 — 0, and elliptic regularity we have 

|(VP((/ • V)g), Vh)\ = |((/ . V)g, (A - VV-) h)\ ^ C\\f\\^, ||Vg||^i \\Ah\y. 

This concludes the proof. □ 

We now return to the proof of Theorem [ 



Proof. We assume there exists a smooth solution u of (|l.ip - (|1.4p on the time inter- 
val [0, T] for some T > 0. We will find appropriate a priori estimates for the norm 
of u (see relation (|2.22l) . below) in terms of the initial data and T. Now a standard 
approximating scheme (e.g. the one constructed in [12 ) will prove global existence 
of solutions. 



Our estimates on the nonlinear term are not optimal. Using 'optimal' estimates here would 
be at the expense of simplicity, and obfuscate the main idea. Further, the 'optimal' estimates are 
still insufficient to prove global existence without a smallness assumption on the initial data. 

''For non-integer values of s, we define the fractional Sobolev norms by interpolation. See for 
instance [H Page 50]. 
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Fix £ > to be small enough so that Proposition 12.21 holds, and {■,■)£ denote 
the equivalent inner product from Proposition [2?2] Then 

(6.4) ^^t\\u\\l^ + {Piiu■^)u),u)^ + {u,Au)^^0. 

where 




By Lemma [6.11 for any cq > 0, we can find a constant C — C{e, cq, f2) > such 
that 

(6.5) \{P{{u-y)u),u)^\ ^ \\u\\^,\\Vu\\H^,A\M\L■^- 

^ CWVutlh^uti^ ^ C\\Vu\f + ^IIA^f. 

We will subsequently fix cq to be the constant c that appears on the right of (|2.7p . 
Using Proposition 12.21 and equations (16.41) . (j6.5|) we obtain 

dtMl. + ^ (llV^^lli. + + a||Vg(«)||i.) ^ C\\Vu\\%, 

where is the constant in (|2.7p . Allowing the constant C = C(£, cq, il) to change 
from line to line, and using the Poincare inequality, we obtain 

dt\\u\\l, + - i\\Vu\\l.. + |||A^||i. + C,\\Vq{u)\\l.^ ^ C\\u\\%, - -Ml,. 

' Cq \ I / = Ci = 

for some constant c\ = ci{e, 17). Thus if at time i = we have 

1 

then for alH > 0, 

\Ht)\\l^ +- f (\\Vu\\l. + |||Au||i. + a||Vg(u)||i.) ds ^ ll^oIlL- 

Co Jo ^ ^ ^ 

Now using the local existence result in [12], and the fact that ||-||jLfi is equivalent 
to the usual norm, we conclude the proof of Theorem 12.61 □ 

7. Two DIMENSIONAL SmALL DIVERGENCE GLOBAL EXISTENCE FOR THE 
EXTENDED NAVIER-StOKES EQUATIONS 

The aim of this section is to prove Theorem l2.7l We recall first the i?Q -orthogonal 
projection onto divergence free vector fields. For uq e H}^{Vl), we define vq^w^ £ 
Hl {n) to be solutions of the PDE's 

{- Awo + V(/> -A7io in f2, ( ~ Awq + VV' = in 17, 

V • Wo = in fl, and < V ■ wq — V ■ uq in 17, 

V = on dfl, I wq —0 on dQ. 

respectively. Note that uq G i?g(17) guarantees the required compatibility condition 
Jj^ V • uo = 0, and so the existence of vo,wo satisfying (|7.ip is well known (see for 
instance [THl §2]). Clearly uo = vq + wq, and orthogonality of vq and wq in 
follows from the identity 

{vo,wo)H^(^n) = (^^0, Vwo) = {vo, -Awq) = {vq, -Vip) = 0. 



22 



GAUTAM IYER, ROBERT L. PEGO, AND ARGHIR ZARNESCU 



Let V be the solution to equation ()l.ip - ()1.4p with initial data vq. As shown 
earlier, V ■ vq = implies that W ■ v — for all time, and consequently w is a 
solution of the 2D Navier-Stokes with initial data vq. Let w — u — v, and observe 

{dtw + P{{w • V)w) + 
+ P{{v ■ V)w) + P{{w ■ V}v) = Aw - Wps{w) in il, 
w ~ on 9ri. 

The strategy to prove Theorem 12. 71 is as follows. First standard existence theory 
for the 2D Navier-Stokes equations implies that for any initial data vq G Hq , with 
V • = 0, we have global existence of a strong solution v. Further, after a long 
time Tg, the solution v becomes small. Now making uq — vq is sufficiently small, 
we can guarantee that w, a solution to (|7.2I) with initial data uq — vq, both exists 
on the time interval [0, Tq], and is small at time Tq. Thus u = v + w is a solution 
to (|l.ip - (|1.4p defined, which is small at time Tq. Now a small data global existence 
result (Theorem 12. 6|) will allow us to continue this solution for all time. 

We begin with a Lemma concerning the existence and smallness of solutions 
to (TT^ . 

Lemma 7.1. Let C 6e a bounded domain, and vq E Hoi^) V-wq — 0. 
Let Uq G HQ{n) be such that PqUq = vq. Then, for any Tq^Sq > there exists a 
(small) constant Wq — Wo{Q, \\vo\\h^ ,To,do) such that if 

Ikollffi ^ Wo 

then there exists a solution of (|7.2|) on the interval [0, Tq] and 

\\w{To)\\h. ^So. 

Momentarily postponing the proof of the lemma, we prove Theorem 12.71 

Proof of Theorem \2.T\ We let Vq be as in Theorem l2.6l and let v be the solution to 
the 2D Navier-Stokes equations with initial data vq = PqUq G Hq- It is well known 
(see for instance [11[T5]) that there exists Tq large enough, so that ||w(To)|Ihi ^ 
iVb- Indeed, from the standard energy identity we can choose To to satisfy 
To{\V,f < ll^olli.. 

By Lemma 17.11 there exists Wq > small enough so that if initially 

(7.3) ll^oll^i ^ Wo 

then the solution w to (|7.2p exists up to time Tq, and further ||w(ro)||^fi ^ 
From (|7.1I) . we know ||wo||^i ^ c||V-wo||^2 (see [15j §2]). Since the norms ||-||^i and 
ll'llj:^! are equivalent, making Uq small enough will guarantee (I7.3p . thus allowing 
to apply Lemma 1 7. II and obtain the existence of w on the interval [0,ro]. Then we 
obtain that (|l.ip - (|1.4p has a solution u = w+v on [0, Tq] and moreover ||w(To)||^i ^ 
||w(To)||^i -I- llwCTo)!!^! < Vq. Applying Theorem 12.61 we can continue the solution 
u on the interval [Tq, oo). □ 

It remains to prove the Lemma. 

Proof of Lemma \7.1\ As with the proof of Theorem 12.61 it suffices to obtain an a 
priori estimate for l|w||^i. Fix e > to be small enough so that Proposition 12.21 
holds. Then 
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Ji 

(7.4) ldt\\w\\l, +'{P{{w ■ V)w),w)^ + {w,Aw)^ 



- {P{{v ■ V)w),w)^ - {P{{w ■ \/)v), w)^ 



where denotes the inner product defined in Proposition 12.21 and ||-||^i the 

induced norm. 

We estimate each term individually. The term J7i is identical to the term that 
appears in the proof of Theorem 12.61 ^-nd thus 

J-i ^ ^ (liv^^ii^. + '-\\Aw\\l. + C,\\Vq\\l.) - (cMl,.^ - ^Ml^^ , 

where q is the unique, mean-zero function such that Vg = (/ — P)w. As before, cq 
is the constant that appears in (j2.7p . and ci = ci{e,fl), C2 = C2(£,r2) are positive 
constants. 

Using C = C(e, f2) > to denote an intermediate constant that can change from 
line to line. Lemma |6. II bounds ^2 and ^Ts by 

IJ2I + IJ3I <C^(||«|lHi||Vii;||^i/.||Aii;||^. + ||u;||^,||Vz;||^i/.||A7«||^.) 

< ^l|Au;||i. + C3 {\\vC. + \\Vv\\l.,.) \\w\\l^ 

for some constant C3 = €3(6, ft). 
Combining our estimates. 



(7.5) U\\w\\l. + - (\\S7w\\l. + j\\Aw\\l. + C,\\S7q\\] 

Z ^ Co \ 4 

Since w is a strong solution to the 2D incompressible Navier-Stokes equations 
with initial data vo € Hq, we have (see for instance [41 p. 78]) that 

/•OO 

snp\\v{t)\\l,+ / \\v{s)\\l,ds<C 
t^o Jo 

for some constant C depending only on 51 and ||t'o|liji . Using this in ()7.5p will prove 
local well-posedness of (|7.2p . Further, for any Tq,So > 0, equation ()7.5p will also 
show that the solution to (j7.2p exists up to time Tq, and ||w(ro)||^fi < Sq, provided 
||wo||j:/i is small enough. □ 

8. Divergence damped equations 

The aim of this section is to prove coercivity of Ba (defined in ()2.25p ). with 
constants independent of a, and 2D global existence with strong enough divergence 
damping (Corollary 12. 8p . 



2 " 



Proposition 8.1. For any a ^ 0, and u G D{Ba) we have 

(8.1) {u,B^u)^^{u,Au)^ + a(\\yQ{u)\\l,+C,\\Qiu)\\l,+e\\AQiu)\\ 
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Proof. By linearity, 

(8.2) (u, Bau)^ = {u, Au)^ + a{u, (I - P)u)^. 

For the second term on the right, 

{u, (I - P)u)^ ^ {u, {I - P)u) + C,{Q{u), QUI - P)u)) + £(Vu, V(/ - P)u). 
The first two terms on the right are equal to ||(/ — P)m||^2 and Ce||Q(M)||^2 respec- 



tively. For the last term, 

\n\ = A(T - P)ii\ 4- / u.. 



Vu, V(/ - P)u) = -{u, A(/ - P)u) + / u,-^ {{I - P)u)^ 



~{u,VW -u) +0 

|V-«|li.- / {V-u)u-,^^\\^-u\\l. 
Jan 



Consequently, 



(8.3) (u, (/ - P)u)^ =\\{I- P)u\\l, + C,\\Qiu)\\l, + e||V • u\\l,, 



and using (|8T2|) . we obtain (|8TT|) . □ 

Before moving to the proof of Corollarv 12. 8[ we digress briefly to remark that 
we can also consider higher order divergence damped operators of the form 

B'^ A - aVV • . 

The results we obtain for (|l.ip - (ll.2p with a zeroth order damping term will also 
apply when we add the second order damping term above. However, while the 
operator B'^ has a stronger (second order) damping term, it is not as easy to deal 
with numerically. The zeroth order damping terms in Ba, on the other hand, can 
easily be implemented numerically, and has a strong enough damping effect to give 
a better existence result (Corollarv 12. 8|) . We now return to prove Corollarv l2.8l 

Proof of Corollarv \2.8[ Since Theorems 12 . 61 and 12 . 71 work verbatim for (j2.24p . there 
exists a time Tq = To(||uo||jyi i independent of a, such that there exists a solution 
u of (|2.24p on the interval [0,To], with ||u(ro)|lHi bounded, independent of a. Let 
Uo be the constant from Theorem l2.7l Observe that (|2.24p implies that V-u satisfies 

dtV ■ u + aV • u = AV ■ u in ft, 
d 

— V-u = foTxedn,t>0, 
ov 

Consequently, 

||V-u(t)||i. <e-(^^+")*||V.^o|li., 

where Ai > is the smallest non-zero eigenvalue of the Laplacian with Neumann 
boundary conditions. Thus there exists ao > 0, such that 

||V.u(To)||^. <C/o, 

for all a > ttQ. Now, by Theorem 12. 71 the solution to (|2.24p also exists and is regular 
on the time interval [Tq, oo). □ 
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9. Existence results under coercive boundary conditions. 

The aim of this section is to show that the extended Stokes operator is coercive 
under the boundary conditions (|2.26p . and prove Proposition 12.91 We begin with 
coercivity. 

Proposition 9.1. // either u and v are in H^{T'^) and periodic, or if u and v are 
in H^{fl) and satisfy the boundary conditions ()2.26p . then 

(9.1) {Au,v) = {u,Av) and {u,Au)^ / |Vu|^. 

Jn 

Proof. In the periodic case, PA — AP. Thus Vps — Q, A — —A, and both 
equalities in (I9.ip follow easily. 

Suppose now u, v satisfy (|2.26p . In view of (|2.ip . we have 

{Au,v) = -{PAu,v) - (VV • u,v) = -{PAPu,v) - (VV • u,v) 

= (VPu, VPv) + {V ■ u,y ■ v) - [ iPv)i^^^^ - (V ■ u)u • 1/ 

Jan I 

Observe that Pv = on dVl^ because because Pv ■ v = Q hy definition of P, and 
Pv • T = by (|2.26p . Thus both the above boundary integrals vanish, giving 

{Au, v) = (VPu, VPv) + (V ■ It, V • v). 

A similar calculation shows 

{Av, u) = {VPu, VPv) + (V • It, V • w) 

proving that A is self adjoint. 

Now because Pu = Pv = on dfl, a direct calculation shows that 

(VPu, VPv) = (V X Pu, V X Pv) = {V X u,V X v). 
Consequently, we see 

{Au, v) ^ {V X u,V X v) + {V ■ u,V ■ v) ^ {Vu, Vv). 
Setting u — V, the second assertion in (|9.ip follows. □ 

Finally, we turn to Proposition 12.91 Before presenting the proof, we remark 
that if we instead impose periodic boundary conditions. Proposition 12.91 and its 
proof (below) go through almost unchanged. The only modification required is 
the justification of the Poincare inequality that will be (implicitly) used in many 
estimates. For this justification, observe that with periodic boundary conditions, 
the mean of solutions to (|1.1|) - (|1.2[) is conserved. Thus, by switching to a moving 
frame, we can assume that the initial data, and hence the solution for all time, are 
mean zero. This will justify the use of the Poincare inequality in the proof. With 
this, we prove Proposition [ 



Proof of Provosition \2.9[ Let v = Pu, and q — Q{u). Since 
P((Vg-V)V<7)-P(V|Vg|V2)-0, 
applying P to (|l.ip gives 

'dtv - PAv + P{{v ■ V){v + Vq) + {Vq ■ V)v) = in Q, 
V = on d^l, 



(9.2) 
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where the boundary condition on v conies from (|2.26p . The point is that energy 
estimates can be used directly to estimate v, since it satisfies expUcit boundary 
conditions. 

Since Pv = v, multiplying (j9.2p by v and integrating yields, 
ld,\\v\\l. + \\S7v\\l^l [ \v\^Aq- I v{{vV)Vq) 

(9.3) < Iw'^vWl. + CMl.wv -uWl^. 

Here we used elliptic regularity to control || V^g|| by ||Ag||, which is valid since ^ — 
QondVt. We also used the (2D) Ladyzhenskaya inequality ||w||^4 ^ C'lklLa || Vw||^2, 
which is valid since w = on (?f2. 

Since V • u is a mean-zero solution of (jl.6p , we know that 



j \\V-u{t)\\l.dt^—\\V-UQ\ 



2 

L2J 



where Ai is the smallest non-zero eigenvalue of the Laplacian with Neumann bound- 
ary conditions. Thus Gronwall's lemma and (|9.3p gives the closed estimate 

(9.4) \\v(t)\\l^ + ds ^ exp (C||V ■ u,\\l,) \\vo\\l.. 

Since Pv = v, regularity of the (standard) Stokes operator tells us that the norms 
||— PAw||^2 and ||w||^2 are equivalent (see for instance i4| Chapter 4]). Multiplying 
(|9.2p by — PAw, integrating by parts, and using (|6.3p to bound the nonlinear term 
in the usual way gives 



dA\Vv\\l. + h^vWl. < C (|kfL2||V^;||i2 + IIVV • ^11^2 + ||V . ut:) \\Vv\ 
Using Gronwall's lemma, equation (|9.4p and (|1.6p . we obtain 

(9.5) W^vml^ + - f\\^<s)\\l. ds ^ KWVvoWl 

c Jo 



2 

L2- 



for some constant K = K{V,, [jV • uol|^2, ||wol|^2). In fact, one can bound K above 

by 

K^Cexp (C (exp (c||V • mWl^) \\vo\\l, + ||V • ^011^2 + ||V ■ uo\\l,^^ 

for some constant C = C{ft). 

Finally, we consider a Galerkian scheme for (ll.ll) ~ (|1.2l) using eigenfunctions of the 
Stokes operator (with no-slip boundary conditions), and gradients of eigenfunctions 
of the Laplacian (with no-flux boundary conditions). It is easy to check that these 
Galerkian approximations satisfy the same energy estimates (|9.4p and ()9.5p . A 
bound for dtu will then follow from (|l.ip . and standard techniques will prove global 
existence. □ 

Appendix A. Failure of coercivity under the standard inner product. 



Most of this section is devoted to the proof that Stokes operator is not positive 
under the standard inner product (Proposition 12. ip . 
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Proof of Proposition \2. 1[ As mentioned earlier, the key idea in the proof is to iden- 
tify the harmonic conjugate of the Stokes pressure as the harmonic extension of the 
vorticity. We begin by working up to this. Since 

Ap, = V • Vp^ = V • (AP - PA) u = 0, 

the Poincare lemma guarantees the existence of Qs such that 

Observe that both ps and qs are harmonic. Indeed, 
(A.2) Ags V X V^qs = V x Vp^ = 0. 

We remark that equations (jA.l[) and (jA.2[) above show that — is the harmonic 
conjugate of Ps- 

To obtain boundary conditions for qg, let t = — be the unit tangent vector 
on dfl. To clarify our sign convention, \i v = (^^), then r {-l-^- Now observe 

^ = Vg,.r = V^g,.^=f^ = z..(AP-PA)« 

= V ■ (Au - VV -u^^v ■ V^V X u = r • Vcj = 

OT 

where, as before, V x m = diU2 — d2Ui is the two dimensional curl, and w = V x u. 
Thus, adding a constant to qs^ we may, without loss of generality assume 

(A. 3) qs ~ io on dfl. 

A direct calculation shows 

u ■ Au — / \Vu\^ + / Vps ■ u — / [Ld^ + lV- u\^ ) ~ / Qs ■ ^, 



In Jfi Jn Jn ^ ' Jn 

where we used the boundary condition u = on 9f2 to integrate by parts. Thus to 
prove Proposition [2T] it is enough to produce a function u, satisfying the required 
boundary conditions, such that 

(A.4) / qsUj^\\uj\\l, + {C+l)\\V -uWl, 

Jn 

We prove the existence of such functions separately. 

Lemma A.l. For any C > 0, there exists u S ^^o(f^) such that (jA.4[) holds. As 
usual, w = V X u, and qs is the solution of the Dirichlet problem 

Aqs — in Q, 

qs ~ itJ on dQ 

The Lemma immediately finishes the proof of Proposition 12.11 □ 



(A.5) 



Proof of Lemma \A.1[ We look for u of the form u = v + Vp with V • t; = in 
and V ■ v — on dfl, where v denotes the outward pointing normal vector on the 
boundary. Then there exists on so that v = V'^tp. The boundary condition 
V ■ I' — becomes ^ = where r denotes the tangential direction on dfl. 
We note that u — and v ■ v — on Oft imply 

dp dip dp ^ dip ^ 

— = -fT = - — , Tp=0, viy ^ — ^0. 
dT dv du dT 
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As ^ = and V' is determined up to a constant we can assume without loss 
of generality that -0 = on dfl and then for a given u the stream function ■0 is 
uniquely determined as the solution of the Dirichlet problem 

Aip — Lu in i7, 
■0 = on dil. 

Then, we have 

qsUjdx— [ qsA->Jjdx= [ Aqgip dx + [ -^q^da— [ -^ijj da 
Jn Jn Jon Jon 

—Qsda^ / —A^da, 
an ov Jg^ dv 

and (|A.4[) becomes 

(A.6) / ^AV'da^ ||A0||^. + (C + l)||Ap|| 

Jon Of 

Summarizing it suffices to find such that (IA.6P holds together with the 
boundary conditions 

(A.7) = 0, f^ = 0, and = on 9^ 

Fix some point xq G dVL and let s i— >■ x{s) be an arclength parametrization of 
the (C^) boundary dil, such that xq = a;(0) and oriented so that the outward unit 
normal j)(s) at x{s) satisfies v{s)^ = x'{s). Then the map {s,r) i— )■ a; = x{s) — ri'{s) 
is and is locally invertible near xq, providing orthogonal coordinates x i— >■ (s, r) £ 
(— e,e) X (0,e) in some neighborhood of xo in ^1. 

We fix p to be of the form p{x) = a(s)/3(r) where a and /3 are in C^((— e,e)) 
and /3(0) = 1, /?'(0) = 0. We will then choose -0 of the form iIj{x) — a'{s)j{r) where 
7 e C;?°((-e,£)) with 7(0) = and 7'(0) = 1. Then (|X7| will hold, and direct 
calculation shows 

/ ^Ai,da= r a'{s)^{j"{0) + K{s))ds 
Jon J^e 

where k(s) = Ar(i;(s)) is the curvature of the boundary. The right-hand side of 
(jA.6p on the other hand, is easily computed to be bounded by C + ^11711^2, with a 
constant C independent of the choice of 7. It is clear that 7 can be chosen to make 
7"(0) arbitrarily large while ||7||^2 remains bounded. Thus (jA.6|l holds for some -ip 
and p. □ 



Finally, to conclude this section we turn to the proof of Corollary 12. 31 Of course 
the proof is immediate from Proposition and we only present it here for com- 
pleteness. 

Proof of Corollarv \2.3\ . Choose mq G C^(f^) to be such that (|2.5p holds, and let u 
be the solution to (|2.10p with initial data mq. By continuity in time, we must have 



u{t) ■ Au{t) < 



n 



for all t in some small interval [0,to]. Thus 9t||M||^2 = ^ j^u ■ Au > on the 
interval (0,to] which immediately completes the proof. □ 
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